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The SU(A^) symmetric antiferromagnetic Heisenberg model with multi-column representations on 
the two-dimensional square lattice is investigated by quantum Monte Carlo simulations. For the 
representation of Young diagram with two columns, we confirm that a valence-bond solid order 
appears as soon as the Neel order disappears at Y = 10 indicating no intermediate phase. In the 
case of the representation with three columns, there is no evidence for both of the Neel and the 
valence-bond solid ordering for N > 15. This is actually consistent with the large-N theory, which 
predicts that the VBS state immediately follows the Neel state, because the expected spontaneous 
order is too weak to be detected. 


Realization of quantum spin liquid in short-range cou¬ 
pling models has been a popular research target in con¬ 
densed matter physics for several decades. One approach 
to obtain a spin-liquid state is to consider a Hamiltonian 
with higher symmetry. Read and Sachdev generalized 
antiferromagnetic Heisenberg into SU(A^) symmetryii^. 
Based on the 1/N expansion they showed that the 
ground state of the model with sufficiently large iV is a 
valence-bond-solid (VBS) breaking the lattice rotational 
or the translational symmetry spontaneously. Recently, 
in terms of the deconfined quantum criticality^Ti^, their 
theory attracts much attention. In particular, the exis¬ 
tence of intermediate state was discussed near the bound¬ 
ary of Neel and VBS^-^. 

Nature of the ground states of the model can vary de¬ 
pending on the representation of SU(7V) algebra. Read 
and Sachdev suggested that for the model with the rep¬ 
resentation of m rows and n columns Young diagram, the 
ground state phase diagram on the N-n plane does not 
strongly depends on the value of m. Within the 1/N ex¬ 
pansion, there are only two types of phases: the small-V 
Neel phase and the large-V VBS phase (see Fig. [TKa)). 
In addition, it was shown that the nature of the VBS 
state can be classified by the remainder of the division 
of n by 4 on the two-dimensional square lattice. For 
n = 1,3 (mod 4), the VBS state is so called colum¬ 
nar VBS where both of translational symmetry and 90° 
lattice rotational symmetry are broken (Fig. HKb)). For 
n = 2 (mod 4), the VBS state is expected to be a ne¬ 
matic VBS with breaking only lattice rotational symme¬ 
try (Fig. me)). In the case of n = 0 (mod 4), there is 
no spontaneous symmetry breaking, which is an analog 
of Haldane state in S' = 1 spin chain. 

Beyond the l/7V-expansion, it was shown that for 
(n, m) = (1,1) the ground state is the SU(V) Neel state 
for TV < 4 while it becomes columnar VBS state for N > 
5 by the quantum Monte Carlo (QMC) calculation^!^. 
Related to this case, the SU {N) J-Q model was proposed 
by Sandvib^. The SU(iV) J-Q model has an additional 
many-body interaction so that the quantum phase tran¬ 
sition between the Neel phase and the columnar VBS 
phase occurs by continuously changing the Hamiltonian 
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FIG. 1. (color online) (a): Schematic phase diagram of the 
SU{N) Heisenberg model on the square lattice with single¬ 
row (m = 1) representations. The phase boundaries for the 
case of n = 2, 3 are determined in the present study, while 
that for n = 1 was from Ref.i^. In the case of n = 3,4, 
we do not see clear evidence of the spontaneous VBS order 
in the vicinity of the phase boundaries for finite-size QMC 
simulations. (b),(c): Schematic picture of the columnar VBS 
(b) and the nematic VBS (c) states. Thick solid lines denote 
larger value of (X )^/3 S/^Sj°‘) while thin solid and dashed 
lines indicate smaller values. 


parameter. In order to vary the Hamiltonian continu¬ 
ously, a continuous-V model was also proposed by Beach 
et al.,—. Because the phase transition might be a realiza¬ 
tion of deconfined quantum criticality^^— , nature of these 
models have attracted much recent interests in condensed 
matter physics^^— . 

For the case of n > 2, there were a few studies con¬ 
cerning phase boundary between the Neel phase and the 
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VBS phase. In QMC calculation up to L = 32 for L x L 
square lattice, no evidence of VBS order was found for 
n = 2, 3,4 with m = 1—. This result appeared to suggest 
an intermediate phase between the Neel phase and the 
VBS phases. However, whether the missing evidence of 
the VBS order for n > 2 is due to an intermediate spin 
liquid phase or due to the extremely small (but finite) or¬ 
der parameter beyond numerical limitation has not been 
clarified up to now. 

In this paper, we investigate ground state of the SU(V) 
Heisenberg model for n = 2 and 3 with to = 1 by using 
QMC simulation. The SU(iV) model we considered is an 
SU(V) symmetric antiferromagnetic Heisenberg model 
on the two-dimensional square lattice with the periodic 
boundary condition. Hamiltonian of the model is given 

by 


«= V V t (1) 

where are generators of SU(iV) algebra, and 

we consider J > 0. On one sublattice A of the lattice, 
the representation of the generators 5'“^ is characterized 
by the Young diagram with a single row (to = I) and 
arbitrary number n of columns, while we use the con¬ 
jugate representation on the other sublattice. Note 
that the conjugate representation satisfies the relation 
We have performed QMC simulation based 
on the loop algorithm. We modified ALPS/LOOPER 
cod o^^d’^ for the present purposei^. We set the inverse 
temperature /3 as /3J = L and investigated the zero 
temperature properties by extrapolating the results to 
L —^ 00. 

In order to see the VBS orders, we define two order 
parameters. The local nematic order parameter is defined 
as 


— Pi,v ( 2 ) 

where (/i = ±x, ±y) is the nearest-neighbor product 
of “magnetic” moments 


N 

Tj _ \ naa naa /o\ 

W.M = 

a=l 

The nematic order parameter characterizes the symme¬ 
try breaking of 90 degrees lattice rotation. ($j) takes a 
finite value for both of the nematic VBS and the colum¬ 
nar VBS states in the thermodynamic limit. We also de¬ 
fine a local complex order parameter characterizing the 
columnar VBS order as 

vk, = (-1)^^ (P,,, - (P,.,^ - P,. _^), (4) 

where jx and jy are integers representing the lattice coor¬ 
dinates of site j. In the columnar VBS phase |(4'j)| ^ 0, 
while |(4'j)| = 0 for the Neel and the nematic VBS 
phases. In order to see the phase transition clearly, 




FIG. 2. (color online) Log-log plot of the two-point correlation 
functions at |Ji| = L/A for the model with n = 2 . (a) The 
correlation function of the Neel order, (b) The correlation 
function of the nematic order. 



FIG. 3. (color online) The binder ratio of the nematic order 
parameter for the model with n = 2 for 16 < L < 128. 


we examine the two-point correlation functions of an 
observable O: Co[R) = (0(0)01 (i?)). For the Neel 
order, we use the correlation of a magnetization 
Cyiag^R) = Sa=i (-R)- We also consider the ne¬ 
matic VBS correlation ONem = 0$(i?) and the columnar 
VBS correlation Cco\{R) = C^{R). 

First, we examine the case of n = 2. In Fig. [5J we 
show the two-point correlations for n = 2at |il| = LjA 
for various sizes L and N. For the case of the Neel or¬ 
der (Fig. [5Ka)), we see that the correlation exponentially 
decays to zero by increasing P for V > 10 while it con¬ 
verges to a nonzero value for N = 9 indicating that the 
Neel state is the ground state for TV < 9 and it is not 
for N > 10. These observations are consistent with the 
previous QMC calculation^. For the nematic order pa¬ 
rameter (Fig. HKb)), the two-point correlations tend to 
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FIG. 4. (color online) Log-log plot of the two-point correlation 
functions at |Ji| = Z//4 for the model with n = 3 . (a) The 
correlation function of the Neel order, (b) The correlation 
function of the columnar VBS order. 


converge to nonzero values for N > 10, although the 
situation at TV = 10 is rather subtle because of larger 
statistical errors comparable with the correlation func¬ 
tion itselfi^. In order to confirm the appearance of the 
nematic VBS order at N = 10, we plot the Binder ratio 
of the nematic order parameter in Fig. [31 The Binder 
ratio for the nematic order parameter is given by 


_ 1 
5Nem = 2 



( 5 ) 


where $ is the sum of local nematic order parameters: 
$ = ffNem is normalized so that gNem = 0 

for the Neel phase, while g^iem = 1 for the nematic (or 
the columnar) VBS phase. The nematic binder ratio at 
iV = 10 develops as the system size is incre ased and, as 
we see in Fig.jSl the order parameter y'($2) slightly devi¬ 
ates upward from the power-law decay, ($2)(L) oc 1/L, 
which should be obeyed asymptotically when the system 
is gapped. It indicates the nematic VBS order at V = 10 
in the thermodynamic limit. We also checked that the 
order parameter 'I' shows no evidence of long range order 
for the case of n = 2. From these observations, we con¬ 
clude that for the case of n = 2 the ground state is the 
Neel state for N < 9, while it is the nematic VBS state 
for N > 10. There is no intermediate phase. 

Next we move to the case of n = 3. We plot two-point 
correlation functions for n = 3 at |i2| = L/4 in Fig. |4l 
For the Neel order, we clearly see from the curvature of 
the curves that the Neel state is the ground state for 
N < 14, and it is not for N > 15. On the other hand, 
we do not see clear difference among different values of 
N in the two-point correlation function of the columnar 
VBS order (see Fig.jlKb)). The behavior of the columnar 
VBS order parameters indicates that the VBS order is 



FIG. 5. (color online) Log-log plot of the size dependence of 
the VBS order parameters for (n, N) = (2, 10) and (n, N) = 
(3,15). The dashed lines represent the fitting curve assuming 
1/L decay. Arrows indicate the correspondence between n = 
2 and n = 3 cases. 


too small to be visible even if it exists for L < 128 finite 
systems with the present statistical errors. Indeed, based 
on the 1/N expansion Read and Sachdev proposed that 
the amplitude of the VBS order parameter becomes expo¬ 
nentially small by increasing N as |('I'j)| ~ ex-p{—NEc) 
with the action of a hedgehog instantoni. The constant 
Ec has been calculated as Ec = cln^ with c = 0.12459 ... 
in the limit V —>■ oo with the spin correlation length ^ 
large but fixed^. 

By using the result of the large N theory^, we try to 
estimate expected amplitude of the VBS order. For the 
columnar VBS order parameter, more precise expression 
for our definition of T'j is given by 

Kd/,)| = ^exp(-7VF;,), (6) 

where a is an unknown constant depending on n/N. For 
the nematic VBS order parameter we also obtain 

|($,)| = ^exp(-VF;,). (7) 

We focus on the expected phase boundary of the colum¬ 
nar VBS phase TV = 15 with n = 3. The spin correlation 
length at this parameter is calculated as ^ ~ 5.2 from a 
fitting of the correlation function of the Neel order. In 
the same ratio of n/N = 0.2, the spin correlation length 
at V = 10 with n = 2 is estimated as ^ ~ 4.7. By substi¬ 
tuting the values of ^ and N into two equations (0) and 
© with Ec — 0.12459In^, we obtain a relation 

|('I'j)Ar=l5,n=3| — 0.67 |($j)Ar=l0,n=2| • (8) 

In Fig. [SI we plot the system size depen dence of the 
columnar VBS order parameter y^(l^'p), where If = 
along with that of the nematic VBS order 

parameter y^These order parameters are expected 
to converge into |('I'j)| and |(d>j)|, respectively, in the 
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thermodynamic limit. For the purpose of better com¬ 
parison, we divide ^ (|4tp) by the factor 0.67 which ap¬ 
peared in Eq. ([8|). In the case of ^ (I'I'P) it decrease as 
expected for the case of no long range order, while 
the L dependence of \J ($2) changes from L~^ around 
L ~ 100 indicating development of a weak long range 
nematic VBS order. From comparison between (|4tp) 
and \/ (d>2), we expect that a signature of the columnar 
VBS order for n = 3, V = 15 becomes visible for the sys¬ 
tems size larger than L ~ 400. Therefore the fact that 
we did not observe any evidence of the long range VBS 
order in the present calculation upto L = 128 does not 
necessarily indicate the presence of intermediate phase 
where both of the Neel and the VBS order disappears. 
Because the QMC calculation for L ~ 400 needs larger 
computational cost than the available resources we can¬ 
not reach a clear answer for the phase boundary in the 
case of n = 3. 

In summary, we have investigated the ground state 
property of an SU(iV) symmetric antiferromagnetic 
Heisenberg model on the two-dimensional square lat¬ 
tice for the representations with the n = 2 and n = 3 
columns Young diagrams. For n = 2, we conclude that 
the Neel state is the ground state for iV < 9 while the 
nematic VBS state becomes the ground state for N > 10. 


Thus, there is no intermediate state between them (see 
Fig. [Ha)). For n = 3, the ground state for V < 14 is the 
Neel state and it disappears for N > 15. Although we ob¬ 
served no evidence of the expected columnar VBS order 
for N > 15, this observation does not exclude the colum¬ 
nar VBS order in this case, because we estimated that 
the signature of the VBS order was invisible for smaller 
sizes L < 400 even if it eventually converges to a finite 
value. Determining the VBS phase boundary for n > 3 
requires further studies. Our analysis indicates that we 
need careful extrapolations of the finite-size data into the 
thermodynamic limit. Naive extrapolations may lead to 
an incorrect characterization of the intermediate region 
even if a weak VBS order is eventually stabilized^S. 
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